Coherent control of steady-state entanglement in driven cavity arrays 
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We show that coherent control of the steady state entanglement between pairs of cavity-atom 
systems in an array of lossy and driven coupled cavities is possible. The cavities are doped with 
atoms and are connected through waveguides, other cavities or fibers depending the implementation. 
We find that the steady state entanglement can be coherently controlled through the tuning of the 
phase difference between the driving fields. It can also be surprisingly high in spite of the pumps 
being classical fields. Furthermore, the maximal of entanglement for any pair is achieved when their 
corresponding direct coupling is much smaller than their individual couplings to the third party. 
This effect is reminiscent of coherent population transfer through a third level between otherwise 
uncoupled atomic levels using classical coherent fields. The above results are analyzed for a range 
of values of the system parameters and different network geometries. 
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Coupled cavity arrays have recently been proposed as 
a new system for realizing schemes for quantum com- 
putation ll| and for simulations of quantum many-body 
systems |2j. More recently driven arrays were considered 
towards the production of steady state entanglement Q 
under realistic dissipation parameters. Also, an anal- 
ogy with Josephson oscillations was shown and the many 
body properties of the driven array have been recently 
studied [|]. 

In this work we examine for the first time the possi- 
bility of achieving coherent control of the steady state 
entanglement. We show explicitly that for a closed ge- 
ometry for a three-cavity setup, the latter being realiz- 
able in a variety of cavity QED technologies including 
photonic crystals, circuit QED and Fabry-Perrot cavi- 
ties connected through fibers @, H, 0, HI, such control is 
possible (see Fig. Q]). Light from the waveguides can di- 
rectly coupled to the photonic modes of the defects cav- 
ities through coherent tunneling. Using the projection 
operator method and following a similar analysis to that 
m Ref. HQ, we can derive an effective master equation 
for polaritons in three cavities. 
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where crj is the polaritonic creation operator for the 
doped cavities, a, the photon field operators of the 
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FIG. 1: (color online). Schematic representation of three- 
doped-defect cavity setup based on a possible implementation 
using photonic crystals (for illustration purposes only): Each 
waveguide carrying the driving classical fields can be replaced 
by fibers in the case of Fabry-Perot cavities or stripline mi- 
croresonators for circuit QED implementations [a, la, 13, lq| • 



central empty cavity and 7, n the dissipation rates of the 
polaritonics and the cavity photons. 
The effective Hamiltonian is 
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FIG. 2: (color online). The concurrence between cavity 2 and 
cavity 3 of Fig. 1 as a function of x\ and y\. This plot is 
similar to the Fig. 2 in Ref. 0]. 

where Ti = J£k/(k 2 + A 2 ) and Aj is the detuning be- 
tween the central-cavity or a fiber coupler Si is the detun- 
ing between the driving field of the central cavity in each 
of the three arms and the polaritonic frequency, assum- 
ing different frequencies for the three cavities, = Gifii, 
Gi is the coupling of the driving field to the central cav- 
ity field in each arm. on is the amplitude of the driving 
field. Note that the photon hopping strengths Ji can be 
complex numbers, however, the phase of Jj does not af- 
fect the entanglement between polaritons, as it can be 
incorporated into the operators of polaritons and thus 
amounts to a local unitary operation. Therefore, we will 
only consider Jj as a real number. 

We can now derive the steady state p ss by requiring 
that p' ss = in Eq. {!]). This is done numerically due to 
the large numbers of coupled equations involved. Next, 
we calculate the entanglement between cavity 2 and cav- 
ity 3. For this calculation, we first need to trace out the 
degrees freedom of cavity 1 and then evaluate the con- 
currence of density matrix of the remaining cavities, 
using 

C(p ss ) = max {0, Ai - A 2 -A 3 - A 4 }, (3) 

where A^'s are, in decreasing order, the nonnegative 
square roots of the moduli of the eigenvalues of p ss p S s 
with 

Pss = {<r v x al)p* a (<rl4), (4) 

and p* s denotes the complex conjugate of p ss - 

The value for the concurrence between cavity 2 and 
cavity 3 (denoting it as C(pj) is plotted against the di- 
mensionless parameters xj = otj(Aj — in)/(JjK), yj = 
Aj/k and Zj = 1 + j/(2Tj). If we assume Ai = A 2 = 



FIG. 3: (color online). The concurrence between cavity 2 
and cavity 3 as a function of </?i and </?3. x\ = 1.67e lvi , 
xz = 1.67e llp3 . When \ip\ — <ps\ = (2k + l)ir (k is an integer), 
the concurrence reaches maximum which is 0.417. The upper 
left figure is the section at 953 = 0. 



A3 = 1.5 x 10 14 Hz and k = 10 13 Hz the maximum con- 
currence can reach 0.402 at z\ — 2:3 = 1.113, z 2 = 114, 
and J/i = J/2 = 2/3 = 15. These correspond to field 
amplitudes ol\ — S3 = 1.215 x 10 3 , and couplings 
d = -G 3 = 1.0 x 10 s Hz, G 2 = & 2 = 0, 7 = 10 a Hz, 
Ji = ,h = 1.0 x 10 12 iiz, J 2 = 3.16 x 10 w Hz. The 
effective dissipation rates appearing in the initial Mas- 
ter equation (Eq. 1) are Ti = T 3 = 4.42 x 10 s Hz and 
T 2 = 4.41 x l0 5 Hz. Theses values are consistent with 
the parameters used in current or near future technolo- 
gies, e.g. with coupled toroidal microcavities quan- 
tum dots [6] and circuit QED ^. In Fig. [2J we plot 
C p for a range of parameters that satisfy x\ — ~x$ and 
J/i = J/3 and find that the behavior to be similar to the 
case of two cavity setup in Ref. [|| , albeit with a much 
higher maximum possible entanglement here. We also 
note here that T 2 << T\ = Ti, and T, cx J 2 , which 
means that the coupling between the two cavities in 
question is much weaker than the coupling between each 
one and the third cavity. Also the state of the polari- 
ton in cavity 1 for the maximum entanglement point is 
Pi = Tr 2 , 3 [p\ - 0.99| ff 0)( 5 0| + 0.01|1-)(1 - | » \g0)(g0\, 
which means that the state for the polariton in cavity 1 
is almost a pure state at ground energy level. Therefore, 
it is almost uncorrelated to the polaritons in cavity 2 and 
cavity 3 which means p « |(/0)((/0| <X> p2,3- Although this 
result initially looks counter-intuitive, it can be explained 
as follows: the maximum entanglement between the two 
parties, i.e. cavities 2 and 3, in a three party system, is 
attained when the state of the third party, i.e. cavity 3, 
nearly factorizes in the combined three party state. The 
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fact that this is happening for strong relative couplings of 
J12 = J\ and J13 = J3 compared to J23 = J2 is reminis- 
cent of the behavior of a coherent process taking place. 
One could dare to observe an analogy here with the case 
of coherently superposing two initially uncoupled ground 
states in a A type quantum system through an excited 
state using two classical fields to mediate the interaction. 

The last observation is further justified by observing 
that C p is larger when x\ = — X3 and y\ — 1/3. Taking 
into account that x\ and X3 are functions of the phases 
of driving fields, these conditions mean that the max- 
imum entanglement is established when the two driv- 
ing fields have opposite phases. In Fig. [3] we plot the 
concurrence against the phases of driving fields. Here, 
z\ = Z3 = 1.01, Z2 — 11. When the phase difference is 
|yi — ^3 = (2fc + 1)71" (k is an integer), we get again a 
maximum of 0.417 and an oscillatory behaviour charac- 
teristic of the expected coherent effect taking place. In 
simple words, when the two fields are completely out of 
phase the entanglement is maximized whereas at phase 
difference 7r/2, the two polaritons are completely disen- 
tangled. In all other cases, the amount of entanglement 
lies somewhere in between. 

In Fig. [31 we remove cavity 1 from the configuration 
of Fig. [1] but with the fiber couplers intact. We then 
analyze the entanglement between cavity 2 and 3. Opti- 
mization of this entanglement gives similar values of the 
parameters like the ones used above (choose z = 1.01); 
however, the concurrence can reach a maximum of 0.47. 
Again the dependence \ip± — 953 1 = (2k + 1)tt (k is an in- 
teger) is apparent (see Fig. [5|. However, if we compare 
the insets in Fig. [3] and Fig. [5] for the cross-sectional 
plots of the concurrence for (p 3 = 0, we see that the plot 
in Fig. [3] has a narrower peak whereas the plot in Fig. 
[5] is broad. This implies that the maximum concurrence 
for configuration in Fig. [4] is substantially more stable 
against variation in the phases (fx and ips than that in 

Fig.ru 

At this juncture, it is worth emphasizing that we now 
have three different configurations for comparisons: (i), 
two cavities with a single driving field in Ref. [|[ ; (ii) , two 
cavities with three driving fields as in Fig. 0] (iii), three 
cavities with three driving fields as in Fig. [1] As shown in 
Fig. 6 when z ranges from 1 to 1.221, the maximum con- 
currence for configuration (ii) decreases rapidly from 0.48 
to 0.285. This rapid decrease means that although con- 
figuration (ii) can reach higher entanglement than config- 
uration (i), yet it is more fragile to the dissipation of the 
environment (parameterized by 7 in z). In comparison, 
the three-cavity setup (which is made by connecting the 
two side empty cavities with each other through a third 
cavity doped with an atom) is more robust against the in- 
crease of dissipation (only when z > 4.03, its maximum 
concurrence drops to be same to that of configuration 
(i)). Therefore, we conclude that cavity 1 in Fig. Q]not 
only coherently mediate between cavities 2 and 3, it also 
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FIG. 4: (color online). Schematic diagram of the two-cavity 
case in which there are three waveguides carrying the three 
respective classical laser fields. Note that each waveguide car- 
rying classical fields can also be replaced by fibers in the case 
of Frabry-Perot cavities or stripline microresonators for cir- 
cuit QED implementation as before. 




FIG. 5: (color online). The concurrence between two cavities 
-Fig. 4- as a function of ^2 and ips,. xi = yi = 0, X2 = 
5e lM ,x 3 = 5e I,fi3 , T 2 = T 3 = 1.316 x 10 8 and Ti = 10 10 . 
When \<p2 — 933 1 = (2k + l)n (k is an integer) , the concurrence 
reaches maximum which is 0.470. The upper left figure is the 
cross section at <p3 = 0. 



stabilized the amount of entanglement betwen the two 
cavities. 

One could try to employ entanglement witnesses to 
detect this entanglement [12|. A witness could be con- 
structed from the density matrix corresponding to the 
maximum value of the concurrence Q and then measure 
it along the corresponding spin directions. In this sys- 
tem to implement the necessary effective spin measure- 
ments we can use the usual atomic state measurement 
techniques employing external laser tuned to the corre- 
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FIG. 6: (color online). The maximum concurrence versus z 
in three configurations: (i), two cavities with a single coupler 
in Ref. |3j; (ii), two cavities with three couplers as in Fig. 
21 (iii), three cavities with three couplers as in Fig. Q] The 
solid/dashed line is for configuration (i) /(ii) when z < 1. 221. 
The dash dot line is for configuration (iii) when z < 4.034. 
The double dot dash line is for configuration (i) and (ii) when 
z > 1.221. The dot line is for all the three configurations 
when z > 4.034. 

sponding polaritonic levels [l3j. In these measurements 
the correlations between the polaritons are transferred 
to emitted photons and can thus be detected by analyz- 
ing the fluorescent photon spectrum. In the following 

we plot the cross-correlation coefficient ^ 2 2 3 f' 3 ^ for 

the three cavity scheme in Fig. 1 as function of the phase 
difference between the driving field 1 and 3 (Fig. [7]). The 
plot is consistent with the concurrence plot in Fig. [3] 

In this work, we showed that steady state entanglement 
in a lossy network of coupled light-matter systems can be 
coherently controlled through the tuning of the phase dif- 
ference between the driving fields. We also find that in a 
closed network of three cavity-atom systems the maximal 
of entanglement for any pair is achieved even when their 
corresponding direct coupling is much smaller than their 
couplings to the third party. This effect is reminiscent of 
coherent effects found in quantum optics such coherent 
population transfer between otherwise uncoupled levels 
through a third level using two classical coherent fields. 
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